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ABSTRACT 

We study a multipole vector-based decomposition of cosmic microwave background 
(CMB) data in order to search for signatures of a multiconnected topology of the 
universe. Using 10 6 simulated maps, we analyse the multipole vector distribution on 
the sky for the lowest order multipoles together with the probability distribution 
function of statistics based on the sum of the dot products of the multipole vectors for 
both the simply-connected flat universe and universes with the topology of a 3-torus. 
The estimated probabilities of obtaining lower values for these statistics as compared 
to the 5-year WMAP data indicate that the observed alignment of the quadrupole and 
octopole is statistically favoured in a 3-torus topology where at least one dimension of 
the fundamental domain is significantly shorter than the diameter of the observable 
universe, as compared to the usual standard simply-connected universe. However, 
none of the obtained results are able to clearly rule out the latter (at more than 97% 
confidence level). Multipole vector statistics do not appear to be very sensitive to 
the signatures of a 3-torus topology if the shorter dimension of the domain becomes 
comparable to the diameter of the observable universe. Unfortunately, the signatures 
are also significantly diluted by the integrated Sachs- Wolfe effect. 

Key words: cosmic microwave background-cosmology: observations 



1 INTRODUCTION 

According to General Relativity, a pseudo-Riemannian man- 
ifold with signature (3,1) is the mathematical model of 
spacetime. The local properties of spacetime geometry are 
described by the Einstein gravitation al field equations. H ow- 
ever, they are only weakly linked ([Roukema et al.l 120071 ) to 
the global spatial geometry of the universe, i.e., its topology. 
We do not have, so far, any theory which could determine 
the topology of our universe, therefore it can be constrained 
only by observations. 

The relea s e of WMA P data ([Bennett et al.1 120031 : 
Hinshaw et al. | I2007L [2009) has stimulated new studies 
of multiconnected (incorrectly dubbed as "non-trivial") 
topology recently, particularly since the detected anoma- 
lies in the observed cosmic microwave background (CMB) 
anisotropy on large angular scales - the putative sup- 
pression of the quadrupole moment, alignment of the 
quadrupole and octopole, and an asymmetry in the sta- 
tistical prop erties of the northern and southern ecliptic 
hemispheres ( de Oliveira-Costa et al. 2004; Copi et al. 2004 
lEriksen et al. 2004; Hansen et al. 2004; Schwar z~eTai1l200- ' 
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- suggest that our universe may possess multiconnected 
topology. 

Amongst all of the multi-connected three-dimensional 
spaces, flat spaces, and the most natural compact mani- 
fold, the 3-torus, have been studied the most extensively 
in a cosmological context. This is motivated by current con- 
straints on the curvature radius of the universe as well as 
computational simplicity. One of the first related analyses 
(|Tegmark, de Oliveira-Costa. &; Hamilton! l2003h suggested 
that a toroidal universe where the smaller dimension was of 
order half the horizon scale may ex plain the above anoma- 
lies, a lthough a more detailed study (|de Oliveira-Costa et al.l 
l2004h did not confirm this hypothesis. Another topology, 
the Poincare dodecah edron for sligh t ly po sitively curved 
space, was proposed bv lLuminet et al.l (|20Q3h . Further stud- 
ies of this topological model revealed some disagreement 
about whether the topology coul d be excluded on the ba- 
sis of the WMAP data or not. lAurich. Lustig &; Steinerl 
(|2005al lbL l2006h : ICaillerie et aD <|2007h claimed that the 
topology can not be excluded, mainly because of insuffi- 
cient accuracy of the data and degradation of the signal 
by the integrated Sachs- Wolfe (ISW) and Doppler contri- 
butions. Moreover, iRoukema etaD (2004) even reported a 
hint of a detection of pairs of matched circles in the CMB 
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maps which fit quite well to t he pr edicti ons for this topol- 
ogy. However, ICornish et all (|2004h and iKev et al.1 (|2007T ) 
ruled it out via the lack of detectio n of statistically signifi- 
cant pairs of matched circles (see also lLew &; Rouke ma 2008; 
iRoukema et~ai1l2008T ). 

Topological studies of CMB maps generally focus on 
two signatures of multi-connectedness: the large scale damp- 
ing of power in the direction of the shorter dimension 
of the doma in, which causes a brea k down of statisti- 
cal isotropy ( Ha nan & Souradeep 2003; iKunz et all 12006. 
120081 : iNiarchou J affe 2006, 200^ 

and the distribution 
of matched pattern s (iLevin et al .11 19981: ICornish et al.lll998l: 
Bond et al.1 l2000allbl: ICornish et al.1 120041: IRoukema et al.1 



2004: lAurich. Lustig fe Steinerl l2005al lbl: IKev et al.1 l2007n 
Levinl (|2QQ2h presents a review of the CMB-related meth- 



ods used in studies of topology. It is important to notice 
that the latter signature is present only when the topologi- 
cal scale is smaller than observable universe, thus methods 
based on it are less powerful than methods based on the 
former signatures of multi-connectedness, which are in prin- 
ciple present also if the topological scale is slightly greater 
than the particle horizon. 

In the former case, one usually uses t he covariance ma- 
trix o f the spherical harmonic coefficients ([Kunz et al . 2006, 
120081 ) of the sky map as the basis of a statistical study. How- 
ever, lAurich et all ()2007l ) used also the so-called multipole 
vectors, initially proposed bv lCopi et al.l ([2004) as a new tool 
for the analysis of CMB maps. In contrast to the coordinate 
dependent coefficients of the spherical harmonic decompo- 
sition, the vectors associated with a given multipole point 
toward the same direction on the sky independently of the 
reference frame employed. Thus, they are useful tools for the 
study of statistical isotropy (irrespective of its topological 
origin, as we consider here, or otherwise). In this paper, we 
study in more detail the application of the multipole vector 
formalism to the detection of signatures of multiconnected 
topology in the context of the observed alignment of the 
quadrupole and octopole moments. Because our studies have 
a rather preliminary character, we will confine ourselves in 
analysis to the simplest compact top ology, i.e., 3-torus t hat 
remains to be ruled out by the data ([Aurich et al.l 120081 ). 

In the following two sections we briefly introduce the 
multipole vectors and related statistics used in the analysis. 
In Section 3] we describe the specific topologies studied. The 
results of the analysis comparing the 5-year WMAP data to 
the studied topologies are presented in Section [5] Finally, 
we summarise the results and draw some conclusions. 



2 MULTIPOLE VECTOR DECOMPOSITION 

The multi pole vector form alism was introduced to CMB 
analysis bv lCopi et al.l ([20041 ), who showed that a given mul- 
tipole moment can be represented in terms of £ unit vectors 
and an overall magnitude, thus making as ex pected 2^ + 1 
figures. As later pointed out by IWeeksl (120041). the formal- 
ism was in fact first discovered long ago bv lMaxwelll(|l89lh . 
Each multipole Tg may therefore be uniquely expressed by £ 
multipole vecto rs v^' 1 ^, . . . , y^'^ and a mag nitude A w . In 
the notation of ICopi et al.l ([2004) , this reads 



Ti(e) = a em Y em (e) = #(v (y) -e) . . . (v (v) -e)+Q , (1) 

where e is the radial unit vector in spherical coordinates. 
Strictly speaking, the multipole vectors are headless, thus 
the sign of each vector can always be absorbed by the scalar 
A w . We will use the convention that all vectors point toward 
the northern hemisphere. 

It is hard to find any convincing physical interpreta- 
tion of the multipole vectors, because they do not indicate 
any specific features of the multipoles. An exception to this 
rule arises when one considers the contribution to a multi- 
pole given by the projection on the last scattering surface 
of a single plane wave. In this case, all the mu ltipole vectors 
point towards the direction of the wavevector (|Lachieze-Revl 
l2004h . Apart from this very specific case, for the quadrupole 
one can notice that the multipole vectors point towards the 
inflection points, such that their cross product points at the 
saddle points. For the higher order multipoles, the vectors 
indicate at vicinity of the inflection points rather than ex- 
tremes. 

To find the multipole vectors from a given set of spheri- 



cal harmonic coefficients, 



one has to solve a set of non- 



linear equations (|Copi et al.1 120041 ). For a quadrupole the 
solution of the equations can be written explicitly (see ap- 
pendix [A}. One should notice that the multipole vectors are 
related non- linearly to the data. For higher order multipoles 
the equations are solved numerically. In our codes we hav e 
implemented the algorithm proposecQ bvlCopi e t al. (2004 
An alternative was suggested by iKatz &; Weeks! (12004 



3 MULTIPOLE VECTORS STATISTICS 

Having computed the multipole vectors, we seek to test a 
given CMB data set with respect to inter-scale correlations 
between di f ferent mul tipoles. In order to d o so, we follow 
ICopi et al.1 (|2QQ4T ) and ISchwarz et all (|2004h and define the 
following set of simple statistics. 

The statistics are based on the dot product, which is 
a natural measure of vector alignment. However, because 
the multipole vectors are defined only up to a sign, one 
needs to consider the absolute value of the dot product. 
Moreover, the multipole vectors do not have their own 
identity, therefore all the dot products were summed for 
a given pair of the multipoles £\ and £2. At the end, to 
normalise the statistics to one, the sum were divided by 
number of the dot products. We consider also the dot prod- 
ucts of the unnormalised cross products of the multipole 



vectors, w 



x ^, k ) as we \\ normalised cross 



products, ,l) = (v 



)/|v ( 



(where 



j / k, j,k = l,...,£ and i = 1,...,£(£- l)/2). Therefore, 
we could use the following types of statistics for any two 
multipoles £\ and £2 (£1 / £2)' 



1 



i x (-2 

EE 

i=i j=i 



^1,0 .^2 ,3) I 



(2) 



The routines of ICopi et al.l (I2004T ) are available at http URL 
http: / /www. phys.cwru.edu/projects/mpvectors~ 
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£ 1 £ 2 (£ 2 -l)/2 

*w = E Iv^^-w^l, (3) 

iWvW i=l j = l 

£ 1 (£ 1 -l)/2 £ 2 

Swv = £ ^w^-v^'l, (4) 

il// wv . . 

1 = 1 J = l 

£ 1 (£ 1 -l)/2£ 2 (£ 2 -l)/2 

Sww = -U £ E l* ( ' 1,0 -* ( ' aJ) l, (5) 

iW ww . . 

1=1 J = l 

where for the last three types of statistics we used also un- 
normalised cross products containing more information than 
the normalised products. Then, the statistics will be denoted 
by iSy W , <$wv an d ^wwj respectively. The four types of the 
statistics are as usual referred to as "vector- vector" , "vector- 
cross", "cross- vector" and "cross-cross" statistics, respec- 
tively. M is number of dot products used for a given statistic 
(M vv « 2 , M vw « 2 (^2 - l)/2, M wv ^1(^1 - l)/2, 
M ww = ^1(^1 - 1)^2(^2 - l)/4, respectively). With such nor- 
malisation the statistics take values in the range [0, 1]. 

Because the statistics are based on the dot product, 
they are rotationally invariant and are sensitive only to the 
relative orientation of the multipole vectors on the sky, and 
not on the absolute orientation of the fundamental domain 
with respect to a given reference frame. 



4 TOPOLOGY OF UNIVERSE 

A detailed description of all possible topologies of three- 
dime nsional m anifolds wit h con stant curvature were given 
bv IWolJ (|l967h and llnoud ((2001). A description and sum- 
mary of the topologies for flat universes, which we will con- 
sider in this paper, and methods of simulat ions of CMB 
maps f or su ch universes were also given by iRiazuelo et al.l 
(|2004aJ ) and IRiazuelo et all (|2004bh . We will briefly describe 
here, following their formalism, basic information concern- 
ing the topology of the 3-torus. 

All flat spaces are obtained as the quotient E 3 /T of the 
three-dimensional Euclidean space E 3 by a group T of sym- 
metries of E 3 that is discrete and fixed-point free. There 
are 18 such spaces among which is the standard Euclidean 
space with trivial topology. Among the 17 remaining ones, 
the ten compact flat spaces are quotients of the 3-torus: six 
are orient able and the rest are non-orient able. The 3-torus 
is the quotient of Euclidean space under the action of three 
linearly independent translations Ti, T2 and T3. The rect- 
angular 3-torus is generated by mutually orthogonal trans- 
lations Ti = (La, 0,0), T 2 = (0,L y ,0) and T 3 = (0,0, L z ) 
where the allowed wave vectors of the cosmological perturba- 
tions k takes the form k = 2ty (n x / L x ,n y / L y ,n z / L z ), where 
n x ,n y ,n z £ Z and L x , L y , L z are dimensions of the funda- 
mental cell, which in this context is called the fundamental 
domain. 

Therefore, there are two effects of the multiconnected 
topology as compared to the Euclidean space: the non- 
isotropic distribution of modes and the discrete power spec- 
trum of the curvature perturbations. Both of these will man- 
ifest themselves in the CMB maps mainly by the ordinary 
Sachs- Wolfe effect, which dominates on large angular scales. 
The two other contributions to the maps, the Doppler and 
integrated Sachs- Wolfe (ISW) effects, will rather dilute the 
signatures of the topology. The former is significant on scales 



smaller than the dimensions of the fundamental domain and 
the latter, though significant on large angular scales, arises 
from the evolution of structure close to the observer. 

To see pronounced signatures of the topology we con- 
sidered domains with at least one of the dimensions smaller 
than the diameter of the observed universe. The distance to 
the horizon is approximately equal to the distance to the 
last scattering surface (LSS) t]lss, therefore the diameter is 
D « 2r/LSS, but depends on the parameters of the cosmo- 
logical model. For the concordance ACDM model, used in 
our work, with the Hubble constant Ho = 70 kms _1 Mpc _1 , 
the matter density parameter Q m = 0.3, the cosmolog- 
ical constant density parameter Qa = 0.7, the diame- 
ter is D = 6.36 Rh, where the Hubble radius Rh = 
c/Hq. We analysed the rectangular prisms with dimen- 
sions L x — L y — 2Rh and L z = 8Rh (hereafter referred 
to as T228) and L x = L y = 8R H and L z 2R H (here- 
after referred to as T882). For such topologies the allowed 
wave- vectors are of the form k = 27r(4n x , 4n y , n z )/L z and 
k = 27v(n x , n y , 4n z )/L z , respectively, where n x ,n y ,n z G Z. 
In Section 15.21 we also show some results for the toruses 
with dimensions L x — L y — 4Rh and L z — 8Rh (hereafter 
referred to as T448) and L x = L y = 6Rh and L z = 8Rh 
(hereafter referred to as T668). The longer dimensions are 
significantly bigger than the diameter of the observable uni- 
verse (i.e., 6.36 Rh) therefore, in practice the topologies are 
indistinguishable fro m the so-called the chi mney and slab 
spaces, respectively (|Adams fe ShapirollioOlh . 



5 RESULTS 

We do not have an analytical expression for either the dis- 
tribution of the multipole vectors on the sky for a multi- 
connected universe or for the probability distribution func- 
tion (PDF) of the statistics. Therefore we based our analysis 
on simulated CMB maps for each of the studied topologies 
with the same cosmological parameters, and in particular 
use 10 5 - for the distribution of the vectors on the sky, and 
10 6 - for the PDF of the statistics. The signatures of a multi- 
connected topology are most clearly visible for the low order 
multipoles, so we confined the analysis only to the multi- 
poles of order t £ [2,8]. For such a range of multipoles we 
have 21 different pairs of multipoles for the multipole vector 
statistics. We did not add noise in the simulations, because 
the WMAP detector noise is completely negligible for the 
low order multipoles. The observed values of the statistics 
are g iven for the ILC WMAP 5 years map ([Hinshaw et al.l 
l2009h . 

5.1 Distribution of the multipole vectors 

In the case of a simply-connected universe we do not ob- 
serve any specific direction in the distribution of the CMB 
anisotropy. As a consequence, the multipole vectors (as well 
the cross products) are uniformly distributed on the sky 
(se<0Fig.|TJ. 

Conversely, in the case of a multi-connected universe, if 

2 because the multipole vectors do not have their own identity, 
in all figures there are shown distributions of one of the vectors 
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Figure 1. The distribution of the multipole vectors (in orthographic projection) in logarithmic scale for the simply-connected universe. 



the dimension of the fundamental domain is comparable to 
the universe's horizon scale, one can see preferred directions 
in the distribution of the CMB anisotropy on the sky. Corre- 
sponding multipole vectors, as well the cross products, also 
demonstrate non- isotropic distributions. The structures we 
observe can be qualitatively explained as follows. For mul- 
tipoles of order £, the biggest contribution comes from the 
modes with the wavenumber k ~ ^/^lss- Significant con- 
tributions will also come from the projection on the sphere 
of those modes with wavenumbers k > ^/?7lss. Because the 
multipole vectors are not related linearly to the spherical 
harmonic coefficients and the map of anisotropy (see ap- 
pendix [A}, it is hard to predict a priori how they are dis- 
tributed on the sky in the case of a multiconnected topology. 
We can, however, expect that the distribution will somehow 
reflect the symmetries of the fundamental domain. 

The distribution of the multipole vectors on the sky for 
the T228 and T882 topologies without contribution from the 
ISW effect are shown in Fig. [2] and Fig. [3] respectively. The 
distribution is related to symmetries of the fundamental cell 
as expected. It is seen particularly for the T228 topology 
and multipoles of order £ = 2 and £ = 4. In the case of the 
quadrupole, the vectors prefer directions pointing towards 
those edges of the rectangular prisms, where the power of 
the maps is concentrated. For the T228 topology these are 
the edges parallel to the XY-plane (in the upper and lower 
sides of the rectangular prism), for the T882 topology, the 
edges perpendicular to the XY-plane. Because the multi- 
pole vectors for the quadrupole point towards the inflection 
points of the CMB maps, we can deduce that the points 



for a given multipole. The distributions of the other vectors are 
the same 



are distributed similarly. The distribution for higher order 
multipoles is less related to symmetries of the fundamental 
cell. It is more azimuthally symmetric and less concentrated 
in the directions preferred by the topology than is found 
for the quadrupole and octopole. It is worth noticing that 
the volume of the fundamental cell for the T228 topology is 
smaller than for the T882 topology. Thus, the non-uniform 
distribution of the vectors is clearer for the former. 

The signatures of the multiconnected topology will be 
diluted by the integrated Sachs- Wolfe (ISW) effect. This 
gives a contribution to the large angular scales. However 
since it arises from the evolution of structures close to the 
observer, it does not contain information about the global 
properties of the universe. The influence of the ISW effect 
on the distribution of the multipole vectors is seen in Fig. 
[4] and Fig. [5] In both cases, the characteristic pattern in 
the distribution of vectors is less pronounced than for maps 
without the ISW effect. 



5.2 The probability distribution function 

The dot products of the vectors for two different multipoles 
are uniformly distributed. However, the distribution of their 
sum does not correspond to the distribution of the sum of 
M independent uniformly distributed random variables, be- 
cause the multipole vectors correspondin g to the same mul- 
tipol e are not statistically independent ( Land &; Magueiiol 
2005). This is seen in Fig. [6] The internal correlations of the 
multipole vectors cause the PDF of the statistic to be nar- 
rower than the distribution of the sum of independent vari- 
ables. Moreover, although at first glance the distributions 
appear Gaussian, they are not so. The kurtosis of the dis- 
tributions is positive, while the skewness is slightly (though 
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Figure 3. The distribution of the multipole vectors (in Mollweide projection) in logarithmic scale for the CMB maps of the multi- 
connected universe with the T882 topology without the ISW effect. 
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Figure 5. The distribution of the multipole vectors (in Mollweide projection) in logarithmic scale for the CMB maps of the multi- 
connected universe with the T882 topology with the ISW effect (cosmological constant density parameter Q/^ = 0.7). 
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significantly) negative for all pairs of multipoles. Deviation 
from Gaussianity is even stronger in the case of universes 
with a multiconnected topology. 

Because the PDF of the multipole vector statistics 
for the multiconnected universes is similar to that for the 
simply-connected universe for higher order multipoles, for 
the T228 and T882 topologies we show only the distributions 
for those lower order multipoles with the biggest deviation 
from the PDF of the simply-connected universe. 

For the T228 topology the alignment of the multipole 
vectors in the direction of the longer side of the prism causes 
a significant shift of the PDF of the S vv statistic toward 
higher values (see Fig. [7}. Because cross products of the 
multipole vectors are orthogonal to the vectors, alignment 
of the multipole vectors results in the distributions of the 
SVw and S wv statistics to be shifted toward zero (see Fig. 
[8] and Fig. [9}. On the other hand, alignment of the cross 
products in the direction of the shorter sides of the prism 
and in the plane perpendicular to the longer side means that 
the distribution of the *S WW statistic is slightly moved toward 
one (see Fig. [TO)) . 

Similar tendencies in the PDF are seen also in case of 
the T882 topology (see Fig. El E21 IE9 and . It can be 
explained by the alignment of the multipole vectors in the 
plane spanned by the longer sides and the cross products 
in the direction of the shorter side. However, because of the 
larger volume of the fundamental domain than for the T228 
topology and weaker alignment of the vectors, the signatures 
of the topology are less significant. Only in the case of the 
5Vw, £ W v and SWw statistics, for the pair of multipoles of 
order £ — 2 and £ = 3 (without contribution from the ISW 
effect), does the distribution demonstrate a bigger deviation 
from the PDF of the simply-connected universe. It has a 
pronounced local maximum around one for the *S WW statistic 
(see Fig. I14[) and show some excess close to zero for the SVw 
and 5 W v statistics (see Fig. [T2l and [T3]) . 

The statistics corresponding to the unnormalised cross 
products take smaller values in comparison to those with 
normalised ones, because of the smaller norm of the unnor- 
malised cross products. It is seen especially well for higher 
order multipoles. 

As in the case of the distribution of the vectors on the 
sky, the ISW effect significantly dilutes the signatures of the 
topology for all kind of statistics. 

For a few pairs of multipoles, especially for (2, 3), (2, 4), 
(3,4) and (3,7), the observed values of the statistics dif- 
fer from the predictions for the simply-connected universe 
substantially. However, not all of them indicate a multicon- 
nected topology. While correlations for the pairs (2,3) and 
(3, 7) prefer one of the topologies under consideration, cor- 
relations for the pairs (2, 4) and (3, 4) are in conflict with 
these models. A quantitative estimation of the goodness of 
fit of the models to the data is presented in Section 15.31 

To study the dependence of the PDF of the statistics on 
the dimensions of the fundamental domain, we estimated the 
PDF also for the 3-toruses T448 and T668. The distribution 
of the statistics for the pair of multipoles (^1,^2) = (2,3) 
without any contribution from the ISW effect is shown in 
Fig. [151 It can be seen that traces of a multiconnected topol- 
ogy are substantially degraded for a universe with larger di- 
mensions, and almost completely disappears already for the 
topology T448. For the simulations including the ISW effect, 



the statistics are even less sensitive to the dimensions of the 
fundamental cell. 

5.3 Probabilities 

To quantify how well a universe with different topologies fits 
the data, we estimated the probabilities of obtaining values 
for the statistics smaller than the data. Since in a multi- 
connected topology the statistics for different pairs of the 
multipoles are correlated with each other (see correlation 
matrices on Fig. [16] and Fig. [T7)) . to correctly estimate how 
well the theoretical models fit we need to consider a joint 
PDF of the statistics for a few pairs of the multipoles. 

From a set of 10 6 simulations we were able to esti- 
mate the joint distribution for up to 4 pairs of the statistics. 
We chose pairs which are the most strongly correlated with 
each other. The correlation matrices, shown in Fig. [16] and 
Fig. 1171 reveal significant correlations between the statis- 
tics where one of the multipole in the pair is quadrupole or 
octopole. The biggest correlations, up to 0.4 for the maps 
without the ISW effect and the T228 topology, are between 
the statistics for the pairs (2,£) and (3, where £ takes 
value in the range £ G [4,8]. They are consequence of the 
alignment of the multipole vectors in the direction of the 
longer sides of the fundamental domains. As was shown in 
Section 15.11 in the case of the quadrupole and octopole the 
alignment is stronger than for higher order multipoles. Let us 
emphasize that these alignments are somehow unrelated to 
the multipole coefficient correlation matrix (a^ m a|/ m /). Ac- 
tually, for parity reasons, coefficients of the second and third 
multipoles are statistically uncorrelated for any torus (i.e., 
(fl2mfl3 TO /) = for all m, m). However, because of the fun- 
damental domain shape, both the quadrupole and octopole 
are led to show up configurations which feel the elongated 
or flattened direction of our torus, and as a consequence, 
their corresponding multipole vector do show some align- 
ment, even though the multipole coefficients themselves are 
uncorrelated. It happens that here, the correlations caused 
by the alignment dominate over correlations between mul- 
tipoles sharing the same parity, like for instance £ = 2 and 
£ = 4, which could be expected to be quite strong. 

We decided to analyse triplets of the statistics for the 
pairs (2,3), (2,£) and (3, <£) (hereafter referred as the £-th 
triplet) as well quartets - for the pairs (2,3), (2,4), (2,5), 
(2,6) (hereafter referred as the quadrupole quartet) and 
(3,4), (3,5), (3,6), (3,7) (hereafter referred as the octopole 
quartet). We chose also a quartet consisting of the pairs 
(2,3), (2,4), (3,4), (3,7) for which the observed values of 
the statistics deviate substantially from the predictions of 
the simply-connected universe. The probability of obtaining 
values at most as extreme as for the data under the null 
hypothesis (i.e., simply-connected universe) should be un- 
usually big for this particular quartet. 

The probabilities for the triplets and quartets of the 
statistics under different hypotheses about the topology are 
presented in Tableland Table [2] respectively. In the case of 
the statistics with cross products of the multipole vectors, 
we present probabilities corresponding to both normalised 
and unnormalised cross products. 

With few exceptions, the estimated probabilities indi- 
cate that universes with a multi-connected topology are a 
better fit to the data, irrespective of whether the ISW effect 
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Figure 6. Distribution of the multipole vectors statistic for the simply-connected universe (black line) and the distribution of the sum 
of M (number of dot products in the corresponding statistic) independent random variables (red line). There are also shown fitted to 
the histograms normal distributions (solid black and red lines, respectively). 
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Figure 7. Distribution of the S vv statistic for the simply-connected (black lines) and the multi-connected universe with the T228 
topology without (red line) and with the ISW effect (blue line) (the cosmological constant density parameter Q\ = 0.7). Vertical line 
indicates value of the statistic for the ILC WMAP 5 years map. 



is taken into account. The T882 topology is preferred by the 
triplets {(2,3), (2,4), (3,4)} and {(2,3), (2,7), (3,7)} as 
well the octopole and {(2,3), (2,4), (3,4), (3,7)} quartets 
for all types of statistics. The preference for the T228 topol- 
ogy by the data depends more strongly on the type of the 
statistic and presence of the ISW effect. 

It is interesting that a few of the triplets of the 
SVv and 5ww statistics, for example {(2,3), (2,4), (3,4)} 
or {(2,3), (2,5), (3,5)}, more strongly favour the multi- 
connected topology for maps with the ISW effect included 
than excluded. This is the case for both of the considered 
topologies. This rather contradicts expectations - the ISW 



effect tends to dilute the signatures of multiconnected topol- 
ogy, therefore the probabilities should be higher for the maps 
without the ISW effect, like for the simply-connected uni- 
verse. It may indicate that the observed relative relations be- 
tween the vectors for these multipoles do not fit well enough 
to the characteristic distribution of the vectors for those 
topologies shown in Section \5. II They only follow the direc- 
tions preferred by the topologies. 

The probabilities for the quadrupole quartet and the 
statistics with the ISW effect and normalised cross products, 
except for the 5 W w statistic, show a stronger preference for 
the T228 topology as opposed to the T882 one. However, 
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Figure 8. Distribution of the SVw and S™ w statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T228 topology without (red line, only for the SVw statistic) and with the ISW effect (blue and light green 
lines, respectively) (the cosmological constant density parameter fl\ = 0.7). Vertical lines indicate values of the statistics for the ILC 
WMAP 5 years map: solid line for S v 
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Figure 9. Distribution of the S wv and S™ v statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T228 topology without (red line, only for the S wv statistic) and with the ISW effect (blue and light green 
lines, respectively) (the cosmological constant density parameter = 0.7). Vertical lines indicate values of the statistics for the ILC 
WMAP 5 years map: solid line for SWv and dashed line for S™ v statistic. 
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Figure 10. Distribution of the S ww and S™ w statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T228 topology without (red line, only for the 5 WW statistic) and with the ISW effect (blue line and light 
green lines, respectively) (the cosmological constant density parameter ^a = 0.7). Vertical lines indicate value of the statistics for the 
ILC WMAP 5 years map: solid line for *S WW and dashed line for S™ w statistic. 
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Figure 11. Distribution of the 5 VV statistic for the simply-connected (black lines) and multi-connected universe with topology T882 
without the ISW effect (red line) and with the ISW effect (the cosmological constant density parameter Q\ = 0.7) (blue line). Vertical 
line indicates value of the statistic for the ILC WMAP 5 years map. 
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Figure 12. Distribution of the S vw and S™ w statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T882 topology without the ISW effect (red line, only for the 5 VW statistic) and with the ISW effect (blue 
and light green lines, respectively) (the cosmological constant density parameter Q\ = 0.7) . Vertical lines indicate values of the statistics 
for the ILC WMAP 5 years map: solid line for SVw and dashed line for S™ w statistic. 
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Figure 13. Distribution of the SWv and S™ v statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T882 topology without the ISW effect (red line, only for the S wv statistic) and with the ISW effect (blue 
and light green lines, respectively) (the cosmological constant density parameter £~2a = 0.7). Vertical lines indicate values of the statistics 
for the ILC WMAP 5 years map: solid line for S wv and dashed line for S™ v statistic. 
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Figure 14. Distribution of the 5 WW and 5 WW statistics for the simply-connected (black and dark green lines, respectively) and multi- 
connected universe with the T882 topology without the ISW effect (red line, only for the 5 WW statistic) and with the ISW effect (blue 
and light green lines, respectively) (the cosmological constant density parameter fl\ =0.7) Vertical lines indicate values of the statistics 
for the ILC WMAP 5 years map: solid line for »5>ww and dashed line for S™ w statistic. 




Figure 15. Distribution for pairs of multipoles (^1,^2) = (2, 3) of the S vv , 5 WW (the left and right figure in the upper row, respectively), 
SVw and S wv (the left and right figure in the lower row, respectively) statistics for the simply-connected (black lines) and multi-connected 
universe with the T228 (red line), T448 (blue line) and T668 (green line) topology without the ISW effect. Vertical line indicates value 
of the statistic for the ILC WMAP 5 years map. 



for the statistics with unnormalised cross products the situ- 
ation is reversed - the data prefer rather the T882 topology. 
On the other hand, for the octopole quartet the statistics 
with normalised cross product prefers the T882 topology 
while with unnormalised cross products, expect for the S^ v 
statistic, - the T228 one is supported. 

In general the probabilities for the statistics are 
higher than for the statistics with the normalised cross prod- 
ucts 5 W . For the multi-connected universe they include a 



few cases with values as extreme as 99% or even 100% for the 
{(2,3), (2,6), (3,6)} triplet. In contrast to the S ww statis- 
tics, they favour a simply-connected universe compared to a 
multi-connected one. 

It is worth emphasizing that, for any of the analysed 
triplets and quartets of the statistics, the computed prob- 
abilities do not clearly rule out the simply-connected uni- 
verse. For the statistics with the normalised cross products 
the {(2,3), (2,4), (3,4), (3,7)} quartet have the biggest 



12 P. Bielewicz, A. Riazuelo 




Figure 16. The correlation matrices of the multipole vectors statistics (S vv , S vw , S wv and 5 WW from the upper to the lower row, 
respectively) for the universe with the T228 topology without the ISW (left column) and with ISW effect (right column). Pairs of 
multipoles are ordered in ascending order of the multipoles i.e., (2, 3), (2, 4), ... , (2, 8), (3, 4), (3, 5), ... , (7, 8). 
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Figure 17. Correlation matrices of the multipole vectors statistics (S vv , SVw, 5 WV and 5 WW from the upper to the lower row, respectively) 
for the universe with the T882 topology without the ISW (left column) and with ISW effect (right column). Pairs of multipoles are ordered 
in ascending order of the multipoles i.e., (2, 3), (2, 4), ... , (2, 8), (3, 4), (3, 5), ... , (7, 8). 
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probability, up to 94%. To a large extent, the correlations of 
the {(2,3), (2,4), (3,4)} triplet are responsible for this be- 
haviour. In this case, similarly large values of the probability 
are seen. For the rest of the triplets and quartets of this type, 
we do not observe any unusually extreme statistical values. 

For the statistics with the unnormalised cross products 
the most extreme probabilities are higher, up to 97% for the 
{(2,3), (2,6), (3,6)} triplet and 96% for the quadrupole 
quartet. 

Before we draw any conclusions from these results, one 
needs to point out that the probabilities shown in Table [1] 
and Table [2] are merely an approximation to the true val- 
ues. Because of the limitation imposed on the number of 
dimensions that we could consider, we did not study the 
joint probability distribution function defined for more than 
4-dimensional space and did not take into account other cor- 
relations, clearly visible in Fig. [16] and Fig. [TT] between the 
statistics. One needs also to remember that the sets of mul- 
tipoles pairs considered were selected on the basis of a pos- 
teriori knowledge concerning the values of the statistics for 
the data. This could also influence our inferences regarding 
the CMB isotropy. 



6 SUMMARY 

We have studied signatures of a multiconnected topology of 
the universe via the distribution of multipole vectors on the 
sky and the PDF of their statistics. We considered topologies 
for a flat universe, where a fundamental domain was a rect- 
angular prism with dimensions L x = L y = 2Rh, L z — 8Rh 
(elongated case) and L x = L y = 8Rh, L z = 2Rh (flat- 
tened case). In both cases, only the shortest dimension of the 
prism is smaller than diameter of the last scattering surface. 
Unsurprisingly, the distribution of the multipole vectors on 
the sky, especially for the quadrupole, reflects symmetries of 
the fundamental cell specific to the topology under consid- 
eration . The multipole vectors are aligned along the longer 
sides of the fundamental cell, but the alignment is more pro- 
nounced for the quadrupole and octopole. As the ISW effect 
is caused by the evolution of structures close to the observer, 
it significantly diminishes the evidence of a multiconnected 
topology, so that the vectors follow instead the distribution 
of the power on the CMB maps. For some of the statistics, 
the joint PDF for a few pairs of the multipoles indicates 
that the data is better fitted by a model of the universe 
with a multiconnected topology. However, one can also find 
statistics for which certain pairs support the alternative con- 
clusions. A more quantitative assessment indicates that the 
data do indeed slightly prefer the multiconnected topology. 
Nevertheless, the most significant values of the probabilities 
for the simply-connected universe, i.e. 94% for the statistics 
with normalised cross products and 97% for the statistics 
with unnormalised products, show that this hypothesis is 
not clearly ruled out. We also found that the multipole vec- 
tors statistics are not very sensitive to the signatures of the 
3-torus topology if the shorter dimension of the fundamen- 
tal domain is comparable to the observable diameter of the 
universe. 

These results have a rather preliminary character. 
There is a need to carefully evaluate the influence of the 
contamination from the Galactic foreground on detection of 



Table 1. The probabilities, Pg, of obtaining values of the statis- 
tics, S vv , Svw, *Swv, *Sww, 5J W , 5J V , 5J W , respectively, smaller 
than the observed values of the statistics for the ILC WMAP 5 
years map. SC denotes the probabilities for the simply-connected 
topology, T228 and T228* - for the CMB maps of the universe 
with the T228 topology without and with the ISW effect, respec- 
tively. Similarly, T882 and T882* denote the probabilities for the 
universe with the T882 topology. 

For pairs of multipoles = {(2,3), (2,4), (3,4)} 





SC 


T228 


T228* 


T882 


T882* 


p s vv 


90% 


92% 


89% 


87% 


88% 


p s vw 


92% 


69% 


82% 


74% 


88% 


p 5wv 


94% 


83% 


90% 


87% 


93% 


P S WW 


94% 


96% 


97% 


84% 


88% 


p ^ w 


73% 


27% 


47% 


65% 


63% 




95% 


99% 


97% 


93% 


94% 




95% 


99% 


98% 


94% 


93% 


?or pairs of multipoles (l 


h,h) = 


{(2,3), (2,5), (3,5); 




SC 


T228 


T228* 


T882 


T882* 


p s vv 


51% 


44% 


17% 


44% 


33% 


p s vw 


69% 


8% 


29% 


36% 


54% 


p s wv 


71% 


12% 


52% 


26% 


63% 


p s ww 


80% 


76% 


79% 


79% 


75% 


p s^ 


74% 


18% 


48% 


49% 


52% 




67% 


68% 


41% 


42% 


52% 


p s- w 


87% 


96% 


98% 


86% 


94% 



For pairs of multipoles (£i,£ 2 ) = {(2,3), (2,6), (3,6)} 





SC 


T228 


T228* 


T882 


T882* 


P 5 VV 


49% 


55% 


26% 


29% 


55% 


P S VW 


86% 


67% 


73% 


63% 


75% 


p s wv 


88% 


20% 


48% 


48% 


69% 


P S W w 


91% 


88% 


88% 


96% 


89% 




85% 


25% 


69% 


64% 


67% 


p s- v 


65% 


56% 


43% 


27% 


47% 


p s- w 


97% 


100% 


96% 


99% 


98% 



For pairs of multipoles (h,h) = {(2,3), (2,7), (3,7)} 





SC 


T228 


T228* 


T882 


T882* 


p s vv 


68% 


40% 


52% 


18% 


40% 


p s vw 


80% 


28% 


39% 


52% 


66% 


p s wv 


82% 


30% 


62% 


62% 


60% 


p s ww 


92% 


92% 


88% 


76% 


78% 


p s- w 


82% 


28% 


44% 


45% 


70% 


P ^ v 


71% 


30% 


37% 


24% 


52% 


P ^w 


92% 


99% 


94% 


97% 


75% 


?or pairs of multipoles (l 


= 


{(2,3), (2, 


8), (3,8); 




SC 


T228 


T228* 


T882 


T882* 


p s vv 


51% 


33% 


34% 


25% 


55% 


p s vw 


77% 


8% 


50% 


51% 


61% 


p s wv 


80% 


27% 


49% 


41% 


64% 


P S WW 


80% 


88% 


76% 


74% 


69% 


p s? w 


81% 


12% 


41% 


37% 


62% 


P S£v 


70% 


66% 


49% 


26% 


49% 


p s- w 


80% 


93% 


93% 


80% 


87% 
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Table 2. The same as in Table ^ but for the quartets of the 
statistics. 

For pairs of multipoles (h , £ 2 ) = {(2, 3), (2,4), (2,5), (2,6)} 





SC 


T228 


T228* 


T882 


T882* 




68% 


oo /o 


*±o /o 


66% 


74% 




88% 


< *± /o 


oo /o 


66% 


86% 




88% 


66% 


81% 


78% 


91% 




89% 


93% 


93% 


80% 


85% 


p s^ 


65% 


29% 


58% 


66% 


56% 




86% 


83% 


81% 


76% 


71% 




96% 


99% 


99% 


85% 


93% 


pairs of multipolee 




= {(3,4) : 


, (3,5), (3,6), (3. 




SC 


T228 


T228* 


T882 


T882* 




88% 


59% 


88% 


79% 


79% 


p s vw 


52% 


48% 


55% 


39% 


51% 


p s wv 


88% 


83% 


78% 


77% 


78% 


p s ww 


88% 


92% 


89% 


67% 


76% 




34% 


15% 


18% 


17% 


29% 




82% 


94% 


87% 


68% 


77% 




80% 


85% 


73% 


70% 


80% 


pairs of multipoles 




= {(2,3). 


, (2,4), (3,4), (3, 




SC 


T228 


T228* 


T882 


T882* 


p s vv 


92% 


87% 


90% 


83% 


84% 


p s vw 


87% 


62% 


76% 


70% 


81% 


p s wv 


92% 


83% 


86% 


86% 


89% 


p s ww 


94% 


95% 


95% 


78% 


83% 


p 5- w 


69% 


22% 


42% 


57% 


62% 




92% 


97% 


95% 


84% 


86% 




91% 


99% 


94% 


89% 


91% 



,7)} 



,7)} 



the topological signatures. It is known that the WMAP ILC 
map used in our studies to some extent is contaminated by 
foreground residuals. The maps for the Q, V and W-bands 
with the masked Galactic plane and corrected foreground 
are more credible for the cosmological analysis. However, 
application of the mask breaks the statistical isotropy of the 
maps. Then, we have to consider all possible orientations of 
the fundamental domain with respect to the mask, which 
makes the analysis computationally much tougher. 

Extension of the studies is also seriously limited by 
the fact that we do not have any analytical expression for 
the PDF of the multipole vectors statistics even for the 
simply-connected universe. So far only the PDF of the mul- 
tipole vectors for the simply- connected universe was given 
([Dennis &; Landll2008l ) . Thus, we had to base our analysis on 
simulations. The correlations between different pairs of the 
statistics present for multi-connected universes introduces 
the necessity of a multivariate analysis of the joint PDF to 
allow a correct estimation of the probabilities. However, us- 
ing 10 6 simulations this could only be achieved up to a four 
dimensional subspace. 
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with corresponding eigenvalues Ai = (X + 3)/6, 

A 2 = A (2) (X - 3)/6, A 3 = -A (2) X/3, where 

X ee v (2 ' 1} • v (2 ' 2) = 3(Ai + A 2 )/|Ai - A 2 | . (A6) 

On the other hand, eigenvector for a given eigenvalue A 
can be expressed in terms of the elements of the Q matrix 

as 

% = ^ [-Qzx (Q xx + Qzz + A) - Q xy Q zy ] , (A7) 

Vy — ~jy [Qzy {Qxx A) Qxy Qzx] 5 (A8) 

V z = jj[{Qxx-\){Qxx + Qzz + \)+Qly\ , (A9) 

where N is normalisation factor chosen such that |V| = 1. 
The eigenvalues are roots of the characteristic polynomial 
A 3 - p A - det Q = 0, where p = Q 2 XX + Q 2 ZZ + Q% + Q 2 xy + 

Qzx Q xxQ zz • 

Thus, using (lA3t . (lA4l) and (lA7|) . (lA8l) . (lA9l) one ob- 
tains direct relation between the spherical harmonics coeffi- 
cients and the multipole vectors 

v (2 ' 1} = Wi V 1 + W 2 V 2 , (A10) 
v (2 ' 2) = WiV 1 -W 2 V 2 , (All) 

where weights Wi y/(l+X)/2 and W 2 y/(l-X)/2 
ensure correct normalisation of the multipole vectors 
|v (2 ' 1} | = |v (2 ' 2) | = 1 and relation (lA6l). 



APPENDIX A: THE MULTIPOLE VECTORS 
FOR QUADRUPOLE 

In case of quadrupole it is possible to show explicitly rela- 
tions between the spherical harmonics coefficients and the 
multipole vectors. Quadrupole may be regarded as 3x3 sym- 
metric traceless matrix Qij such that 

2 

T 2 (e) - J2 fl 2^ m (e) = A (2) [(v (2 ' 1} • e)(v (2 ' 2) • e) + 



- (2,1) 



>(2,2) 



Q ij X X 



(Al) 



Considering Cartesian coordinates x % — (z,x,y) with con- 
1 the Qij matrix takes the form 



straint x 2 + y 2 + z 2 



Qij — \ -x n 



( 



\ 



2 alf 



\/2a\ 



«21 



— y/2 a 2 ™ 



V2a lm 



a 20 
' V3 



a 2 i 
al n 2 



As po inted out bv lLand &: Magueiiol (2005) and lCopi et aU 
(2006) eigenvectors of Q are related to the multipole vectors 
by 



V 1 



OC 
OC 



-(2,1) . -(2,2) 



-(2,1) 



-(2,2) 



at-3 -(2,1) w -(2,2) 



(A3) 
(A4) 
(A5) 



.(A2) 



